Anthropic versus cosmological solutions to the coincidence problem 
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In this paper we investigate possible solutions to the coincidence problem in flat phantom dark 
energy models with a constant dark energy equation of state and quintessence models with a linear 
scalar field potential. These models are representative of a broader class of cosmological scenarios 
in which the universe has a finite lifetime. We show that, in the absence of anthropic constraints, 
including a prior probability for the models inversely proportional to the total lifetime of the uni- 
verse excludes models very close to the ACDM model. This relates a cosmological solution to the 
coincidence problem with a dynamical dark energy component having an equation of state param- 
eter not too close to —I at the present time. We further show, that anthropic constraints, if they 
are sufficiently stringent, may solve the coincidence problem without the need for dynamical dark 
energy. 



I. INTRODUCTION 

More than a decade has elapsed since the first obser- 
vational evidence in favor of cosmic acceleration [l|, 
Since then, a growing body of independent data has con- 
firmed that our universe is indeed accelerating [3,4]. Cur- 
rent observational constraints on the equation of state 
parameter Q suggest that a cosmological constant, A, 
might be playing the role of dark energy. However, the 
value of A required to explain the observed cosmic accel- 
eration is off by ~ 120 orders of magnitude from the stan- 
dard quantum field theory prediction. Another problem 
associated with this model is that we seem to be living in 
a very special epoch where cosmological acceleration has 
just started. This is known as the coincidence problem. 

Dynamical dark energy models alleviate some of the 
problems associated with the cosmological constant. The 
so-called tracker quintessence [1-0] and k-essence mod- 
els consider a scalar field which tracks the be- 
haviour of the dominant component. Also, in interacting 
quintessence models [lfj| - ll3| the dark energy field inter- 
acts non-minimally with the matter fields and suitable in- 
teraction terms may lead to an attractor scaling solution 
in which the ratio between dark energy and matter den- 
sities remains constant. These models have in common 
the feature that, except at very early times, the dynam- 
ics of the scalar field responsible for the dark energy is 
roughly independent of initial conditions. However, they 
provide no convincent cosmological explanation for the 
coincidence of our observing time with the onset of cos- 
mological acceleration. Allowing for several accelerating 
periods has also been considered as a possible solution to 
the coincidence problem [f4l - fl6j . This way, the present 
accelerating epoch would be only one of many. However, 
the major drawback is that there is no evidence that the 
universe has undergone an accelerating period in the past 
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(other than primordial inflation). 

In [l7l ] an upper bound on A was obtained, taking into 
account that larger values would not allow the forma- 
tion of collapsed structures necessary to accomodate life. 
This led to the so called multiverse scenario where the 
cosmological constant is a random variable (one value 
for each of the universes in the multiverse ensemble) . It 
has been argued that anthropic constraints, in such sce- 
nario, select values of A very close to the observed one 
[rsL fig ] . In [2(| an estimate was made for the temporal 
distribution of terrestrial-planet-bound observers which 
has been used as an argument in favor of an anthropic 
solution to the coincidence problem (see also [2l|, [22]). 
However, due to the lack of sufficient understanding of 
the conditions required for intelligent life to appear, an- 
thropic constraints are usually associated with very large 
uncertainties. 



Cosmological models in which matter and dark energy 
densities are of the same order of magnitude for a signif- 
icant fraction of the universe's lifetime avoid the coinci- 
dence problem because in such models the present epoch 
is no longer special. In this paper we consider phantom 
dark energy models with a constant dark energy equation 
of state [U H3] and quintessence models with an effective 
linear scalar field potential potential [25l - t28| . These mod- 
els are representative of a broader class of cosmological 
scenarios in which the universe has a finite lifetime, thus 
allowing for a cosmological solution to the coincidence 
problem if the universe's lifetime is not much larger than 
the age of the universe today. Anthropic constraints, 
parametrized by a time cut-off t c above which observers 
cannot exist, are also taken into account but we do not 
attempt to estimate t c . The main aim of the paper is 
then to determine the values of t c and w for which a 
solution to the coincidence problem is cosmological, an- 
thropic or both. Throughout the paper we shall assume 
the metric signature [+, — , — ] and use units in which 
c = h=87rG= 1. 
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II. DARK ENERGY MODELS 

Consider dark energy models described by a real scalar 
field 4> minimally coupled to gravity with action 
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where £ m and £((/), X) are, respectively, matter and dark 
energy Lagrangians, X — (f>^(f>'^/2 and a comma rep- 
resents a partial derivative. If <fi^ is timelike then the 
energy-momentum tensor of the real scalar field can be 
written in a perfect fluid form 



T» v = (p + p)u^u v - pg^ 



(2) 



with 
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p = 2X£. x -£, p = £(X,<f>), (3) 



In Eq. @, is the 4- velocity field describing the motion 
of the fluid, while p and p are its proper energy density 
and pressure, respectively. The equation of state param- 
eter w is equal to 



P 

w = — 
P 



£ 



2X£, X - £ 
and the sound speed squared is given by 
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£ x + 2X£ XX 
The energy-momentum tensor of the matter field is 

IT = PmVV , (6) 

where v M is the 4- velocity field of the matter field and p m 
is its proper energy density. Its proper pressure, p m , is 
equal to zero so that both the equation of state param- 
eter and the sound speed vanish (w rn — p m /p m = and 

cL = o). 

In a flat homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) universe the line element can 
be written as 



ds 2 



dt 2 - a 2 (t) (dx 2 + dy 2 + dz 2 ) 



(7) 



where a is the scale factor, t is the physical time and 
x, y and z are comoving spatial coordinates. Einstein's 
equations then imply that 



- = --(p(l+3w)+p m ) 
a o 



(8) 



where w = pj p is the equation of state parameter of the 
dark energy and a dot represents a derivative with re- 
spect to physical time. Energy-momentum conservation 
for both matter and dark energy components leads to 



Pro — —3Hp m , 
= -3H(l + w)p, 



(9a) 
(9b) 



where H = a/a is the Hubble parameter. Eqs. (l§a|) and 
(pb|) imply that p m = p m0 a~ 3 while p = p a~ 3 ' 1+ ™^ for 
constant w. The subscript '0' refers to the present time 
which we shall define by the conditions f2,„o = 0.27 and 
do = 1 throughout the paper. Here Q m — p m /Pc and 
p c = 'AH 2 is the critical density (in a flat universe p c is 
equal to the total density pt = p m + p)- 



A. Model I 

The simplest homogeneous and isotropic phantom dark 
energy model has a constant equation of state parameter, 
w < — 1. In [23[ it has been shown what would be the 
fate of such universe. The energy density of phantom 
energy increases with time, pushing the universe apart 
with increasing strength as time goes by. Eventually, the 
universe would reach the point where any strongly grav- 
itationally bounded objects would be ripped off - such 
end of the universe has been dubbed a Big Rip. Mathe- 
matically, this corresponds to a singularity in the FRW 
metric scale factor a(t), which becomes infinite in a finite 
interval of time. The universe's lifetime, is given by [23l | 
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The simplest realization of the above model considers a 
scalar field Lagrangian with a negative kinetic term given 
by 



£ = -X-V(<t>), 



(11) 



Dynamically such models are very similar to quintessence 
models except that, in the phantom case, the scalar field 
climbs up the potential instead of rolling down. The 
sound speed squared is equal to unity and 



-X - V 

-x + v ' 



(12) 



is smaller than —1 as long as V > X > 0. By requiring 
that w < —1 is a constant, one obtains 



V = Vna 



-3(w+l) 



A + B In 



(1 + Vl + ra 3w Y 



(13) 



(14) 



where B = y/Z\l + w\/Zw, A = B In (l + vT+7) with 
r = f2 m o/^eOj so that 4>o = <Xa = 1) = 0. 

Eqs. (|T3"|) and (|14p can be combined to yield the po- 
tential as an explicit function of the field 



V = Vn 



2(1 + VT+OeH 



v 2(l + VT+T)+r(l-e*) / 



(15) 
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It is possible to show that with this choice of potential 
the equation of state of dark energy approaches w inde- 
pendently of the initial value of cj). The same also holds 
in the case of the constant w quintessence models inves- 
tigated in [H-Hi] (in fact Eqs. (13]), flU]) and (15]) also 
apply in that case). The major drawback of this choice 
of potential is that it requires a great deal of fine tuning 
at the transition between the matter and dark energy 
dominated eras to prevent w to change rapidly there (see 
31, 32] for a detailed discussion). For the sake of our 
analysis, it will be sufficient to take this model as rep- 
resentative of phantom models that lead to a future Big 
Rip singularity (sec 33] for a more detailed discussion). 



B. Model II 

Consider a quintessence real scalar field cf> described by 
the Lagrangian 



C = X - V(cf>) , 



(16) 



where the potential is assumed to be linear with V(4>) = 
Vq + V\(j) (Vq and V\ are real constants). The assumption 
of a linear potential is always valid for any realistic model 
of this type for a certain time interval around the present 
day. 

We shall again assume that <f> = today so that Vq is 
the scalar field potential energy at the present time. The 
scalar field equation of motion is given by 



3H(j> = -Vx 



(17) 



We shall also assume that the kinetic energy of the field cj) 
has no memory of initial conditions (28| . This appears to 
be a reasonable assumption given that, in this model, the 
memory of initial conditions is rapidly erased. It has been 
shown (see 26, 28] for a more detailed description) that 
this model has a decelerating matter dominated phase 
followed by a transition to an accelerating dark energy 
phase. As the field continues to roll down the poten- 
tial, eventually V(4>) becomes negative which triggers the 
rapid collapse of the universe into a Big Crunch. This sce- 
nario depends mainly on the linearity of the scalar field 
potential for — Vq < V(cf>) < Vq irrespective of the specific 
form of V{4>) outside this range. For a given Hq, there 
are only two degrees of freedom, Vq and V± , which corre- 
spond to different choices of the parameters f2 rn o an d w 
(we shall reduce the number of degrees of freedom to one 
by fixing fl m Q = 0.27 which essentially fixes Vq). We may 
then use the freedom to vary V\ (and consequently the 
value of the equation of state parameter of the dark en- 
ergy at the present time wq) to numerically compute dif- 
ferent possibilities for the total lifetime of the universe t u , 
using a standard ODE solver. In this model the weighted 
average value of w in the time interval [0,to]j 



Wq 



J Q a ° da (1 - tt m ) u 
J a ° da (1 - Q m ) 



(18) 




FIG. 1: The value of Ht as a function of the cosmological time 
t for four representative cases: w = — 1 (solid line), w — —0.98 
(dotted line), model I with w = —1.1 (dot-dashed line), and 
model II with wo = —0.98 (dashed line). The coincidental 
state, indicated in the figure by the horizontal stripe, was 
defined by Ht — 1 ± 0.05. This defines a coincidental epoch 
indicated by the vertical stripe. 



C. The age of the universe 

The age of the universe (in units of H^ 1 ) as a function 
of a is given by 



Ht = H(a) 



rfa 



aiT(a) 



H(a) f a 



da 



H Jq y/n m0 a, 1 + (1 - fi m0 )i 



-3w-l 



(19) 



satisfies 1.6(1 + wq) ~ 1 + wq [26 



where t(a = 1) = tQ an d w was assumed to be con- 
stant in the last equality of Eq. (T5]). The value of Ht 
evolves from being very close to 2/3 deep into the mat- 
ter era to 2/(3(u> + 1)) at late times if — 1 < w < 0. If 
w < — 1 then Ht — > oo as a — > oo (this happens in a 
finite timescale if w < — 1). If dark energy is modelled 
by a quintessence field with a linear potential, then Ht 
can no longer be computed using the last equality in Eq. 
(IT51) and must be evaluated numerically. In this case, Ht 
increases from 2 /3 at the start of the dark energy dom- 
inated era only to start decreasing at a later time. The 
negative potential energy density of the scalar field then 
triggers the collapse of the universe with Ht — > — oo as 
a — > in a finite time (note that a < during the collaps- 
ing phase) . The fact that the latest observational data is 
consistent with i?o^o ~ 1 is another way of describing the 
coincidence problem. This is illustrated in Fig. [1] which 
shows the value of Ht as a function of the cosmological 
time t for four representative cases: w = — 1 (solid line), 
w = —0.98 (dotted line), model I with w = —1.1 (dot- 
dashed line) and model II with wq — w(a = 1) = —0.98 
(dashed line). Here we have taken the conservative es- 
timate Ht — 1 ± 0.05 in the definition of coincidental 
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FIG. 2: The value of to/t u as a function of |l+too| for model I 
(solid line) and model II (dashed line). The linear approxima- 
tion (dotted lines) given by Eq. (I21|l holds for |1 + Wo | < 10~ 
in the case of model I and |l + wo| < 10 -2 in the case of model 
II. 



FIG. 3: Posterior probability distributions P (wo\0) for 
model I (dashed line) and model II (dot-dashed line) assum- 
ing that a — 0.05 and no anthropic constraints. The results 
are reasonably well described by the analytical approximation 
given by Eq. (|24p (solid line), specially in the case of model 
I. 



state (horizontal stripe in Fig. Q]). This also defines a 
coincidental epoch, At, indicated by the vertical stripe 
(note that At is almost independent of wq for wo ~ — 1). 
Fig. [1] shows that we seem to be living at a special cos- 
mological time, to ~ Hq 1 , very close to the start of the 
dark energy dominated era. If w = — 1 (solid line) then 
t u = oo and consequently Ai/i u = 0. The same applies 
if to > —1 (dotted line). On the other hand, in models 
I and II the coincidence problem may be strongly alle- 
viated. In these models, the coincidence may be much 
less dramatic if the universe ends not too far into to the 
future, either in a Big Rip (w < —1, dot-dashed line) 
or in a Big Crunch (quintessence with linear potential, 
dashed-line) . The coincidental epoch At, corresponding 
to the coincidental state Ht = 1 ± 0.05, may therefore be 
a significant fraction of the total lifetime of the universe. 

III. RESULTS 




FIG. 4: The same as Fig. [3] but with a = 0.005. The analyt- 
ical result given by Eq. (|24[1 becomes a better approximation 
for smaller values of a. 



A. Cosmological constraints 



t u is not too small then 



According to the Bayes' theorem 



P(w \O) 



P (Q\w Q ) P pr i or (w ) 
P(0) 



(20) 



where P pr ior(wo) is the prior probability distribution of 
wo, P (O\wo) is the probability of the data given the pa- 
rameter wq and P (wq\0) is the posterior probability dis- 
tribution of wq (in model I to = mjo)- The term P (O) is 
simply a normalizing constant. In our work, the prior 
probability in Eq. (|20[) is assumed to be proportional to 
the fraction of the total lifetime of the universe corre- 
sponding to a coincidental state where Ht = 1 ± 0.05. If 



t u ~Ct /\l+w \ 



(21) 



is a good approximation in the case of both models de- 
scribed in the previous section. Here C > is a model de- 
pendent constant of order unity (C ~ 0.79 and C ~ 0.13 
in the case of models I and II, respectively). Hence, 



P, 



prior 



(io ) oc 



At 



-|l + itfo|. 



(22) 



where we used the approximation given by Eq. (I21[) . In 
Fig. [21 we plot the value of t /t u as a function of |1 +Wo\ 
for model I (solid line) and model II (dashed line). The 
linear approximation (dotted lines) given by Eq. (|21[) 
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FIG. 5: Posterior probability distributions P (wo\0) for 
model I with t c = oo (solid line), t c /C = 50io (dot-dashed 
line) and t c /C = 150io (dashed line), assuming that a = 0.05. 




FIG. 6: Posterior probability distributions P (wo\0) for 
model II with t c = oo (solid line), t c /C = 50io (dot-dashed 
line) and t c /C = 150to (dashed line), assuming that a = 0.05. 



holds for |1 + ioo| < 10 _1 in the case of model I and 
|1 +n>o| ^ 10~ 2 in the case of model II. For simplicity 
we shall assume that in model I P (0\wq) is given by a 
gaussian distribution with 



P(0\wq) oc exp 



2a 2 



(23) 



if wo < —1 and zero otherwise. In model II we shall 
assume Eq. ([2"B"]) to be valid for wq > — 1 with P(O\w ) = 
if w < -1. Using Eqs. ([20 ]) . (|22 ]i and ([23 ]) . one obtains 
the posterior probability distribution 



P(w \O) 



\l+w \ 



exp 



(1 + wof 
2a 2 



(24) 



which satisfies the normalization condition 
f:^P(w \O)dw = 1 or f™P{w \O)dw = 1 in 
the case of models I or II, respectively. The peak of 
the probability distribution is at |1 + iuq| = cr. Note 
that the constants of proportionality in Eq. (1221) . which 
are model dependent and vary with the definition of 
coincidental state, are absorbed in the normalization 
of Eq. ([24]). In both models, using Eq. ([24]), the 
probability 



P(\l+w \ < A) = 1 - exp 



A 2 \ A 2 
-2^ ~2^' < 25 > 



where the approximation is valid for A/ a < 1 . 

Figs. [3] and H] show the posterior probability distribu- 
tion as function of |1 + wo\/a for models I (dashed lines) 
and II (dot-dashed lines) in the absence of anthropic con- 
straints, assuming that a = 0.05 and a = 0.005, re- 
spectively. Figs. [3J and H also show the analytic esti- 
mate of the posterior probability obtained using Eq. ([24]) 
(solid lines). The quantitative differences registered are 



expected because for small values of t u the linear rela- 
tion in Eq. ([2"T]) . t^ 1 oc |1 + Wq\, is no longer valid (see 
Fig. [5]). The posterior probability distribution P (wq\0) 
differs less from Eq. (l2~4l in the case model I than in 
the case of model II since Eq. (|21l) holds for a wider 
range of 1 1 + Wq \ in the phantom model than in the linear 
quintessence model. Lowering a has the effect of decreas- 
ing the statistical importance of the values of 1 1 + wo \ for 
which Eq. (|21[) is not a good approximation and conse- 
quently both models have posterior probability distribu- 
tions that fit best Eq. ([M]) with cr = 0.005 than with 
cr = 0.05, as can be seen by comparing Figs. [3]and[3] In 
particular, in Fig. [4] the posterior probability distribu- 
tion of model I (dashed line) almost coincides with the 
analytical approximation given by Eq. (I24[) (solid line). 

Eqs. ([2~T1) and ([2"2"]) show that values of wo very close to 
— 1 lead to very large values of t u and, consequently, to 
very small values of the fraction of time where Ht ~ 1. 
The effect of the prior is therefore to strongly disfavor 
wo = —1, even when observations that lean towards a 
cosmological constant are considered. This is clearly seen 
in Figs. [3] and H] which show that the peak of the distri- 
bution is shifted away from wo = —1. This result links 
a cosmological solution to the coincidence problem to a 
significant deviation from the ACDM model. Although 
models with a constant w/ -1 can be constructed, they 
do require a large amount of fine tuning [H [H Hi and 
consequently any deviation from w = — 1 should be inter- 
preted as a hint for a dynamical dark energy component. 



B. Anthropic constraints 

Anthropic constraints may be added to the analysis by 
taking into account that there may be time windows for 
the existence of observers. The corresponding probability 
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FIG. 7: The value of /, defined by Eq. as function 

of t c /C for models I and II assuming that A = 0.05 (solid 
and dotted lines, respectively) and A = 0.005 (dashed and 
dot-dashed lines, respectively). 



distribution function is poorly known (see however [201 ]) 
and consequently in this paper we shall simply assume 
that there is a cut-off at a time t c above which observers 
cannot exist. This way, assuming that the linear approx- 
imation given by Eq. (|2ip holds and C ~ 1, the prior 
probability distribution becomes constant for t u > t c , 
or equivalently |1 + wo\ < to/t c . Hence, the posterior 
probability distribution P {wq\0) no longer vanishes for 
Wq = — 1. On the other hand, for t u < t c , or equivalently 
|1 + w o| > to/t c , the cut-off has no effect on the shape 
of the prior probability given by Eq. (l22j) . As a result 
Eq. (|2~5j) remains valid as long as t c ^> to/ A but it is no 
longer a good approximation for t c < to/ A. Unlike the 
t c = oo case the constant C can now affect the results. 
However, under the transformation t c t c /C, the above 
discussion remains valid. 

The effect of the anthropic cut-off t c on the results 
of models I and II is shown in Figs. [5] and H2 respec- 
tively, which plot the posterior probability distributions 
P(w \O) with t c = oo (solid line), t c /C = 50t (dot- 
dashed line) and t c /C — 150io (dashed line), for a = 
0.05. One sees that below some value of jl+wol ~ Cto/t c , 
the posteriors aquire the gaussian shape of P(0\wo), con- 
sequence of the constant prior. If t c is sufficiently large 
(dashed lines) then the posterior probability distribu- 
tion remains essentially unaltered if |1 + wq\ > Cto/t c , 



becoming nearly constant for lower values of |1 + wo|. 
On the other hand, for smaller values of t c (dot-dashed 
lines) the change of P(wq\0) for |1 + wo\ < Ct /t c with 
respect to the t c = oo case has a significant impact 
on the normalization of the probability distribution for 
|l + tuo| > Cto/t c . 
Fig. [7J shows the value of 



P(|l + ico| < A) t 



(26) 



P(\l + w \ < A) tc=co 
as function of t c /C for models I and II assuming that 
A = 0.05 (solid and dotted lines, respectively) and A = 
0.005 (dashed and dot-dashed lines, respectively). For 
large enough t c the ratio / approaches unity while for 
small values of t c it increases significantly showing that 
the anthropic constraints may have a large impact on 
the analysis. For small enough t c the prior probability 
distribution becomes independent of |1 + wq\ which is 
responsible for the plateaus in Fig. [JJ The differences 
between the results of models I and II for the same A, 
in particular for small t c , result from the breakdown of 
the linear approximation given by Eq. (1211) at larger 
(smaller) values of 1 1 + Wq \ in the case of model I (model 
II). This is also the reason why, in Figs. [5] and [6] for the 
same value of t c /C, the change of P(wo\0), with respect 
to the t c = oo case, is more noticeable in the case of 
model II than in the case of model I. 



IV. CONCLUSION 

In this paper we investigated cosmological as well as 
anthropic solutions to the coincidence problem in the 
context of universes ending in either a Big Crunch or a 
Big Rip, considering two representative families of cosmo- 
logical models. In the absence of anthropic constraints, 
values of Wq very close to —1 are ruled out when a prior 
probability inversely proportional to the total lifetime of 
the universe is included in the statistical analysis. When 
anthropic constraints are taken into account, the analysis 
no longer predict a null probability density for wo = — 1. 
The probability distribution is a function of the assumed 
time cut-off t c for the existence of observers which is un- 
certain. Depending on the value of Wo and t c a solution 
to the coincidence problem can be either anthropic, cos- 
mological or both. 
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